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A time-dependent electric field gives rise to a stationary 
non-equilibrium current I^ 2 ' around a mesoscopic metal ring 
threaded by a magnetic flux. We show that this current, 
which is proportional to the intensity of the field, is closely 
related to the exchange part of the interaction contribution to 
the equilibrium persistent current, and that the corresponding 
non-linear conductivity directly measures the weak localiza- 
tion correction to the polarization. We explicitly calculate the 
disorder average of J 1 - 2 ' in the diffusive regime as function of 
the frequency of the electric field and the static flux piercing 
the ring, and suggest an experiment to test our theory. 

PACS numbers: 73.50.Bk, 72.10.Bg, 72.15.Rn 



I. INTRODUCTION 

Electron-electron interactions in disordered mesoscopic 
metals are not very well understood. The usual per- 
turbative machinery of many-body theory is not always 
applicable in these systems, because often the intricate 
interplay between interactions and disorder in phase co- 
herent systems cannot be described by means of simple 
perturbation theory. The persistent current Q in a meso- 
scopic diffusive metal ring threaded by a magnetic flux 
belongs to this category. Seven years after the seminal 
experiment by Levy et al. g], there seems to be general 
agreement that electron-electron interactions are essen- 
tial to explain the surprisingly large magnitude of the ex- 
perimentally measured average persistent current in an 
array of 10 7 Cu-rings. Note that the experiment by Levy 
et al. has recently been independently confirmed || . 

In this work we shall study electron-electron interac- 
tions in mesoscopic metal rings by means of a somewhat 
unconventional approach, which is based on the connec- 
tion between electron-electron interactions on the one 
hand, and non-linear response to an external electromag- 
netic field on the other hand. In the context of persistent 
currents this connection has recently been pointed out 
by Kravtsov and Yudson [0, who considered the time- 
independent part of the non-equilibrium current propor- 
tional to the intensity of an external longitudinal electric 
field E{t) = Re[E(uj)e lu % 



I {2) = Re \a {2 \u,(j)) \ \E{u)\ 



(1) 



u> of the external electric field, as well as of the static 
flux cj> piercing the ring. The fact that non-linear re- 
sponse and interactions are closely related becomes obvi- 
ous in a path-integral approach. Indeed, it is well known 
H that the Coulomb interaction between electrons can 
be obtained by integrating the exponential of the cou- 
pled Maxwell-matter action over the fluctuating quan- 
tum electric and magnetic fields. However, in a path 
integral approach we may also perform the integrations 
in a different order. Thus, an alternative method to ob- 
tain the equilibrium current is to calculate first the non- 
equilibrium current for a given realization of the elec- 
tromagnetic fields, and then performing an average over 
these fields. The effective action for this averaging pro- 
cedure is obtained by integrating first over the electronic 
degrees of freedom, keeping the electromagnetic fields 
fixed. In this work we shall show that this procedure 
leads to new perspectives in the role of electron-electron 
interactions for persistent currents, which can be tested 
experimentally by means of non-linear transport exper- 
iments. We shall also re-examine and correct an ear- 
lier calculation of the non-equilibrium current (Q) due to 
Kravtsov and Yudson B. 

The equivalence between Coulomb interactions and 
fluctuating electromagnetic fields has been employed pre- 
viously by Altshuler, Aronov and Khmelnitsky || in their 
calculation of the dephasing rate due to electron-electron 
interactions in disordered metals. See also Ref. Q] for a 
recent study of interaction effects in mesoscopic conduc- 
tors with the help of this approach. 

The plan of the rest of this paper is as follows: In Sec.|| 
we shall use well-known functional techniques to show 
how the above averaging procedure can be carried out 
in practice, and elucidate the precise connection between 
non-linear response and electron-electron interactions. In 
particular, we show that the non-equilibrium current (|l|) 
is closely related to the Fock contribution to the equilib- 
rium persistent current, and that the non- linear conduc- 
tivity cr (2) (w,0) in Eq.(|) can be simply obtained from 
the weak locali zati on correction to the polarization of the 
system. In Sec. Ill we shall explicitly evaluate the average 

_(2) 

non-equilibrium current I as function of u and 4>, and 
compare our result with Ref. 0. (Here and below the 
over-bar denotes averaging over the disorder.) In Sec. IV 
we discuss possibilities to test our theory experimentally, 
and conclude in Sec.[v] with a brief summary. 



The so defined non- linear conductivity <j( 2 )(cj,0) (see 
Eqs.(E^) and (Mq) below) is a function of the frequency 
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II. FROM INTERACTIONS TO NON-LINEAR 
RESPONSE 

In this section we shall map the problem of calculat- 
ing the equilibrium persistent current of interacting elec- 
trons onto an effective non-equilibrium problem in imag- 
inary time. This is achieved by means of a Hubbard- 
Stratonovich transformation, a well known technique in 
the theory of strongly correlated electrons |§J| . 



A. Definition of the problem 

We consider a system consisting of electrons with 
charge — e and mass m which are confined to a thin ring 
with circumference L and cross section L\, with L± L. 
The electrons interact with two-body Coulomb forces and 
move in a static random potential U (r) . Identifying the 
position along the circumference of the ring with the x- 
coordinate, the equilibrium current / around the ring can 
be written as 



I = — I drj(r) , 



(2) 



where the integral is over the volume V = LL\ of the 
ring, and the particle current density j(r) (in x-direction) 
can be expressed in terms of the exact imaginary time 
Green's function G(r, r', r — r') as 



j(r) = lim lim J x>x ,G{r, r', r - r' - 0+) 

r' — >r v' r 

oo 

= T V e iC » 0+ lim J X)X ,G(r,r',iw n ) , (3) 

* — * r'^r 

n— — oo 

with the differential operator 

Jx.x' = - — :{d x -d xl )-\ , a=— — . (4) 

Zmi m L <po 

Here 4>o = 27rc/e is the flux quantum, and we have intro- 
duced the imaginary frequency Fourier transform of the 
Green's function, 



G{v,v',t-t')=T 



71— — OO 



G(r,r',iQ n ) , (5) 



where ui n = 2n{n + ^)T. For simplicity, we shall work 
with spinless electrons and use units where Ti and the 
Boltzmann-constant are set equal to unity. This amounts 
to measuring temperatures T and frequencies to in units 
of energy. The Green's function at constant chemical 
potential fi can be represented as a functional integral 
over Grassmann fields ip and ^ in the usual way fe|, 



G(r,r',T-r') = - 



where S{ip} = S {ip} + S'int{V'}, with 

kk' 

SintW = ipf, X] flP-qPq ■ 



(7) 

(8) 



Here p q = J^ fc iplipk+q, and the inverse non-interacting 
Green's function matrix in the momentum-frequency ba- 
sis for a given realization of the disorder potential is 



[G 1 ]fc/ C ' — 6, 



kk' 



(k + a) 2 

luJ n X h /-i 

Am 



— 5 nn i Uk- 



(9) 



where a is a vector potential directed along the circum- 
ference of the ring (which in our convention is identified 
with the x-direction) , with magnitude a = |a| = x"^- 
For simplicity we have introduced collective labels k — 
[k, iQ n ] and q = [q, iuj-n] for wave- vector and Matsub- 
ara frequencies, where ui n — 2n{n + i)T is a fermionic 
frequency, and uj n = 2itnT is a bosonic one. The Grass- 
mann variables ipk are the Fourier components of the 
field V(r»T), i.e. ^(r,r) = V" 1 / 2 J2k e*( k - r - £> » r >Vfc- The 
Fourier transforms C/ q of the disorder potential and / q 
of the Coulomb potential are normalized such that both 
have units of energy [fill, i.e. 



C/ q = ~ J dve-^U{v) , 



(10) 
(11) 



The disorder potential U (r) is assumed to have zero av- 
erage and Gaussian white noise correlations, so that 



UqU, 



q' 



(12) 



where the parameter 7 is a measure for the strength of 
the disorder. Within lowest order Born approximation 
we may identify 7 = A/(2wt), where A is the average 
level spacing at the Fermi energy, and r is the elastic 
lifetime. 

The evaluation of the above expression for the current 
would require the solution of the many-body problem in 
the presence of disorder, an impossible task. Perturba- 
tive expansions can be performed in powers of the disor- 
der potential t/ q and in powers of the Coulomb interac- 
tion / q . This double expansion is rather subtle. To ob- 
tain sensible results which correctly take into account the 
physics of diffusion and screening, infinitely many powers 
of Uq and / q have to be summed. In order to make this 
expansion more transparent and to see the connection 
with non-linear response, we shall now map this problem 
onto an equivalent problem where the two-body interac- 
tion is replaced by a time-dependent auxiliary field. 



(0) 
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B. Hubbard- Stratonovich transformation and 
equivalent non-equilibrium problem 

The two-body part Si nt {^p} of our effective action 
can be decoupled by means of the following Hubbard- 
Stratonovich transformation || 



(13) 



Applying this transformation to the denominator and nu- 
merator of Eq.(^), and integrating over the Grassmann 
field, the exact current density of the many-body system 
can be written as 



/ X>{$}exp 


-7 E,/* 1 *-^ 


/£>{•!>} exp 


-lE./q 1 *-^ 





-Seff{*} 



j(r,r,{$}) 



/P{$}e- S cff{<i»} 
(j(r,r,{$})) Sol , . 



The effective action S'egj^} is given by 



ScffW^E/q^-A-Trm 



1-G V 



(14) 



(15) 



where the matrix elements of the infinite matrix V are 
given by [V] k k> = Vfe-fe' = iT^h-h'- The quantity 
j(r, t, {$}) is the non-equilibrium current density for a 
frozen configuration of the Hubbard-Stratonovich field, 
i.e. 



j(r,r,{$}) = Um 4,*'S(r, r', r, r + 0+) . 

r' — >r 

Here Q satisfies the partial differential equation 
(-«V r + 



(16) 



-9r 



2m 



+ Ii-U(t)-V(t,t) 



£(r,r', t,t') 
(17) 



= 5(r-r')5*(r-r') , 
where the time-dependent potential V(r,r) is defined by 

y(r,T)=£y(* r -" - r V 9 , v fl = iT3 9 , (is) 



and <5*(t) = J , ^] I] e" lk '° T is the antiperiodic imaginary 
time 5-function. Note that the potential V(r, r) is a pe- 
riodic function of r, i.e. V(r,r+ 1/T) = V(r, r). On the 
other hand, the fermionic Green' s function Q has to sat- 
isfy antiperiodic boundary conditions in each imaginary 
time variable pT| , 

Q{v,v',t + 1/T,t') = -g(r,r',T,T') 

= g(r,r',r,r'+l/T). (19) 

The above transformation is exact, and allows us to 
clarify the precise connection between interactions and 



non-linear response Q. Note that Eq.(p"7|) defines the 
imaginary time non-equilibrium Green's function of non- 
interacting fermions subject to an external imaginary 
time potential V(r, t). Of course, for a comparison with 
experiments, which measure the non-linear response to 
external fields, we need to know the real time dynamics. 
While within linear response the well-known fluctuation- 
dissipation theorem tells us how to obtain the real time 
response by simple analytic continuation from an imag- 
inary time formalism, in the case of non-linear response 
the situation is more complicated. Nevertheless, even 
then the analytic continuation from the imaginary time 
response to real times is possible, provided the time- 
dependence of the external potential can be analytically 
continued, and the potential is adiabatically switched on 
]l2| . This point, which apparently is not widely appreci- 
ated in the literature, has already been discussed in the 
classic textbook by Kadanoff and Baym . 



C. How functional averaging reproduces the 
equilibrium current 

For a calculation of the equilibrium persistent current 
to first order in the RPA (random phase approximation) 
screened interaction, it is sufficient to expand the effec- 
tive action (|l^) and the non-equilibrium current-density 
j(r, t, {$}) defined in Eq.(|l6|) to second order in the 
Hubbard-Stratonovich field. The resulting Gaussian in- 
tegrations can then be performed exactly. In this approx- 
imation the effective action is given by 



T 



E[<W/q '+110 (<?,<?')] 



where 



N (q)=5 n0 N (q)=Tj2[Go: 



k-\-q,k i 



no(g, ?') = <wn (q, q', iu n ) 

= -T2_Jfio]k+q,k' +q' [Gp]k'k 



(20) 



(21) 



(22) 



kk' 



Physically iVo(q) is the spatial Fourier component of the 
density, and Ho(q, q', iio n ) is the non- interacting polar- 
ization p0| for a given realization of the disorder po- 
tential Uyr). To expand the non-equilibrium current- 
density, it is convenient to consider the Fourier compo- 
nents, 



(r,r,{$}) = jE ci(q '"" r: 



Oq 



(23) 



Note that with this normalization the equilibrium current 
defined in Eq.(0) is simply given by 
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1 - —(jq=o)s a{{ , 



(24) 



where q = means q = and Lo n = 0. For j q we obtain 
the following expansion in powers of V q — iT$ 9 , 



^4 0) +4 1) +4 2) 



where 



ft 

9' 

4 2 ) = ^^ 2 )(,,g', g ")^V- 



(25) 

(26) 
(27) 
(28) 



The linear response function K^ x \q, q') can be identified 
with the non-interacting correlation function between 
density and current-density, 



i^W) = T]T 

kk' 



k x + a + q x /2 



[Go] k+q,k'+q' [Go]k',k ■ 

(29) 



The quadratic response function is 



K^\q,q\q")=T ]T 



k x + a + q x /2 



[Go]k+q.k 



'+q' 



X [Go] k' ,k" [Go] k"~q",k 



(30) 



Graphical representations of jq 1 ^ and j q 2 ^ are shown in 
Figs.|l](a) and (b). It is instructive to see how func- 
tional averaging of these expressions with the effective 
action (^0|) yields the well-known Jl3|,0 interaction cor- 
rections to the equilibrium current to first order in the 
RPA interaction. Of course, the equilibrium current is 
more easily obtained from the derivative of the thermo- 
dynamic potential with respect to the static flux 0,H , 
but the following calculation clarifies the close connection 
between non-linear response and electron-electron inter- 
actions [Q. To perform the Gaussian integration, it is 
convenient to first eliminate the linear term in Eq.(|20|) 
by redefining the $-field such that its Gaussian aver- 
age vanishes. This is achieved with the help of the 
shift-transformation $ 9 = $ g - iT~ x J2 q > f q ^' A No(q'), or 
equivalently for V q — iT$ q , 



V q = V q + Y,f^ A N {q') 



(31) 



Here f™? A is the inverse of the infinite matrix with el- 

J qq 

ements Sqq'f^ 1 + U (q,q'). It follows that within the 
Gaussian approximation 







(V q V- q ,) Sea = ~Tf^ A 



(32) 
(33) 



Substituting Eq.(|l]) into Eqs.@ and (||), and averag- 
ing over the T^-field in Gaussian approximation, it is now 
easy to show 



(4 X) K« ^EE^t^o/r^fe) . ( 34 ) 

q' qi 

<4 2) >*o ff = EE if(2) fe«'-3")/^ I A /^o(gi)iVofe) 

q'q" 9192 
, \— * r , i ,is pRPA 



q'q" 



(35) 



Graphically Eq. (|34|) and the first term in Eq.(p5|) can 
be represented by the Hartree diagrams shown in Figj^, 
while the second term in Eq.(|35|) is represented by the 
Fock diagram in Fig.[| To see that for q = Eqs.(|34j) and 
( |35| ) reduce to the well-known [jl3 14 first order (in the 
RPA interaction) corrections to the equilibrium current, 
we use the following exact identity 



]TGo(k,kV, 

k' 



A' „ 



m 

L<t> d 
2tt d(b 



-G (k',k",^„) 



Go(k,k",iw r . 



(36) 



where [G ] kk > = 5 nn >G (k, k', iu> n ). Eq.(|36|) can be easily 
proven by taking the derivative of both sides of Eq. ( |l7j ) 
(with V(r, r) set equal to zero) with respect to a. With 
the help of this identity we obtain for the functional av- 
erage of the linear response current (p4|) 



-J-(J0 /Seff 



•iE4T^[^o(-q)iVo(q')] 

qq' 



(37) 



where f^i A = /qo P q'o ^ s ^ ne static RPA interaction. Sim- 
ilarly, for q — Eq.(pq) can be written in the form 



-e,.( 2 )> 

-t(Jq )s ett = -~}^ 



RPA 

qq' 



JV (-q)JVb(q') 



(38) 



where 



qq' 



/-RPA 



^ - f&Mn J^ U o(<l, q', tw„) (39) 



is the Fock contribution to the equilibrium persistent cur- 
rent JPH|]. The sum of Eq.@ and the first term in 
Eq.(|38[) yield the (non self-consistent) Hartree contribu- 
tion to the equilibrium current, 



^ = -|^E/ q T^o(-q)iVo(- 



(40) 



We have argued elsewhere |T^] that the neglect of self- 
consistency in Eq.(40) does not properly take into ac- 
count the subtle interplay between disorder and inter- 
actions, so that the correct order of magnitude of the 
Hartree current can only be obtained by means of a self- 
consistent calculation. 
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III. QUADRATIC RESPONSE TO AN 
EXTERNAL ELECTRIC FIELD 

A. Derivation of the non-equilibrium current from 
the Fock correction to the equilibrium current 

Our rather unconventional derivation of the interac- 
tion correction to the equilibrium current makes the con- 
nection between electron-electron interactions and non- 
linear response manifest. In fact, from our derivation it is 
clear that the Fock contribution (^) to the equilibrium 
current is closely related to the non-equilibrium current 
given m Eq.(§)' Physically our Hubbard-Stratonovich 
field $ can be identified with the scalar potential of elec- 
tromagnetism, which is generated self-consistently by the 
motion of the electrons ||,|| . Therefore the negative gra- 
dient of our auxiliary potential V(r, r) is the effective 
force acting on the electrons, which in turn can be asso- 
ciated with an internal electric field E(r, r), 



- eE(r,r) = -W(r,r) 



(41) 



Defining E(r,r) = Y, q e l(q ' r_w " T) E 9 and using Eq.Q 
we have eE„ = iqV„, or 



V q = -ie 



(42) 



where q = q/|q|. From Eq.(B3) we thus conclude 



RPA 
qq' 



|q||q 



-((q.E g )(q'-E_ g ,)>S eff , (43) 



so that the Fock contribution (|39|) to the equilibrium cur- 
rent can be written as 

I F =^a( 2 )(q,q',z^ l )((q-E 9 )(q'-E_ 9 0)s off , (44) 



where the non-linear conductivity (q, q', iu> n ) is given 
by 

er (2) (q,q',iwr0 = ~, M ,, ^-n (q, q', iuJ n ) • (45) 
2 |q||q'| o<p 

Note that functional averaging restores translational in- 
variance in time, so that the average in Eq. ( p!4| ) is propor- 
tional to S nn i . This equation shows that the Fock current 
can be viewed as the sum of functionally averaged non- 
equilibrium currents, generated in second order in the 
internal electric fields associated with the motion of the 
electrons. Clearly, the corresponding non-linear response 
function o^> is a system property that should be inde- 
pendent of the origin of the electric fields. In particular, 
if we add an external electric field, Eq.(|44|) is still valid 
provided we identify E with the total electric field. Thus, 
after performing in Eq.(45) the usual analytic continua- 
tion, ioj n — * u> + i0 + , we conclude that the non- linear 



conductivity <r^ (u>, <fi) defined in Eq.([l]) can be identi- 
fied with 



cr (2) (w,0)= lip a {2, {q,q',u + iO + ) 



(2) 



(46) 



Recall that a^ 2 \uj , <j)) describes the time-independent 
part of the non-equilibrium current that is proportional 
to the intensity of a time- dependent, spatially uniform 
electric field. 

Two remarks are in order. The first concerns the an- 
alytic continuation of the imaginary frequency response 
to real frequencies. Because Eq.([45|) relates a non-linear 
response function to the flux-derivative of a linear re- 
sponse function (the polarization), we may rely on the 
fluctuation-dissipation theorem to relate real and imag- 
inary time response. Of course, Eq.( j45|) c an also be ob- 
tained with the help of the identity (|36|) from the ex- 
pression ([SO]) for the general non-linear response func- 
tion K^ 2 \q,q' ,q"). It is not difficult to see that quite 
generally the correct real frequency response can be ob- 
tained by replacing iuj n — > u> + i0 + , iu> n > — > lo' + i0 + , and 
iuj n n — » lo" + i0 + for all frequencies (l2) . 

Secondly, we would like to emphasize that so far we 
have not averaged over the disorder potential, i.e. all 
equations given above are valid for an arbitrary realiza- 
tion of U(r). In the following subsection we shall start 
from Eq. ( flq ) to perform the disorder average of the non- 
linear conductivity. This has the advantage that we only 
need to average a product of two Green's functions. Note 
that in the work (J] the average of the non-linear con- 
ductivity has been calculated directly from Eqs.(|28|) and 
( ^p| ) . Because the latter expression involves a product of 
three Green's function, it is in this case more difficult to 
identify the dominant disorder diagrams. One difficulty, 
which apparently has not been noticed in Ref . Q , lies in 
the fact that at long wavelengths and small frequencies 
density-vertices are renormalized by singular vertex cor- 
rections involving so-called diffuson propagators jL4 16 
This might partially explain the discrepancies between 
the result of Kravtsov and Yudson j|] and our result dis- 
cussed below. 



B. The relation between non-linear conductivity, 
polarization, and linear conductivity 

Using the fact that Ilofq, a' , iuj n ) = <5q q 'IIo(q, iui n ), 
we see from Eqs.([44|) and (|45| ) that the disorder aver- 
aged static non-equilibrium current that is generated at 
quadratic order in a space- and time-dependent longitu- 
dinal electric field with Fourier components E(q, lo) is 
given by 



/ (2) =^ 2 >(q,^)|E(q, W )| 



with 



7?( 2 ) 



(q,w,0) = §^2^rn (q,w) 



(47) 



(48) 
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We thus need to know the flux-dependent part of the 
disorder averaged polarization of the ring. For frequen- 
cies \ui\ ;$ A and wave-vectors |q| $ 2n/£ (where I — vpr 
is the elastic mean free path) non-perturbative methods 
are necessary to calculate this quantity |17|. Here we are 
interested in the high frequency regime \u>\ ^ A, where 
we may use the impurity diagram technique Jig] , How- 



for 



u> 



t ;$ 1 and \q\£ ;$ 1 the direct diagrammatic 
calculation of ^IIo(q, to) is not so easy, because there ex- 
ists non-trivial cancellations between vertex corrections 
to the density vertices |l6|,[l4]]. Physically, these correc- 
tions arise from the diffusive motion of the electrons in 
the disordered metal. To take these corrections into ac- 
count without having to perform complicated manipula- 
tions, we use the exact relation between irreducible polar- 
ization and longitudinal conductivity er(q, u) |nj , which 
in our normalization fllCIl reads 



- q 2 V_ 
n(q,w) = i -er(q, w) . 

u> e 

From Eq.(Q) we thus obtain 

er (2) (q,w,0) = £_!i_^CT(q,w) 

2 (— iu>) o<p 



(49) 



(50) 



At finite q, the dynamic conductivity has a diffusion pole 
[pToft . In fact, according to Ref. |^0| in the limit of small 
wave- vectors and frequencies 



cr(q,w) 



V Hq' 



■a(u>) , 



(51) 



where <t(w) = ef(0,w), and the frequency-dependent dif- 
fusion coefficient is related to the dynamic conduc- 
tivity via pQ| 



V{uj) _ a(oj) 
"AT ~ 



00 



(52) 



Here X>o is the classical diffusion coefficient, which is re- 
lated to the Drude conductivity do via the Einstein rela- 
tion "Dq = (AV) _1 e 2 uo- We thus conclude 



a™ 



-v— 

2 86 



2?(w)q 2 — iuj 



(53) 



Diagrammatically, the diffusion pole in Eqs.(|5lj) and (|52^ 
implicitly takes the so-called diffuson diagrams into ac- 
count ]l8| . On the other hand, the weak- localization 
corrections described by the Cooperon diagrams have to 
be included explicitly in the calculation a(uj) and T>(lu). 
These diagrams are responsible for the dominant depen- 
dence on the magnetic flux. 



latter is determined by the famous weak localization cor- 
rection arising from coherent backscattering []l8| , 







cr(w) = — 



ttV 86^ P (Q + 2a) 2 



1 



(54) 



where the prime means that the sum is restricted to 
|Q| ^,2tt/£, and L v is the dephasing length [Q. Scaling 
out the Thouless energy E c = T>q/L 2 and setting now 
q = in Eq.(||), we obtain from Eq.(|J) 



cEcjeL/Egf 
(-iuj)(f) 



9{u,4>) > 



(55) 



where u) = to / E c , and the dimcnsionlcss function g{oj, < 
is given by 

I ,iU-^-V' — 

9[UJ ' P) 2tt 86^ p (Q + 2a)2 - lLU + V /Ll ■ 



(56) 

For a thin ring with L± ^ £ the Q-summation in Eq. ( |5^ ) 
is one-dimensional, and can be carried out exactly, with 
the result 

2e~ w sm(47T0/0o)[l-e- 2W ] 



<p) = 



W [1 - 2e~ w cos(47T0/0o) + e~ 2W ] 



(57) 



Here W — W (L/L v ) 2 — iw, where the root has to be 
taken such that ReW > 0. For \W\ <C 1 this reduces to 



4>) 



4sin(47r0/0o) 



[4 sin 2 (2n6/6 ) - iQ + (L/L v ) 2 ] 



(58) 



Note that by definition L v > L in a mesoscopic sys- 
tem, so that for \uj\ <C E c the parameter |H^| is smaller 
than unity. On the other hand, for |o;| ^ E c the pref- 
actor e~~ w in Eq.([57|) reduces to exp[— ^/|w|/2], so that 
the non-linear conductivity becomes exponentially small. 
We disagree in this point with Kravtsov and Yudson [Q, 
who found that the non-linear conductivity in the regime 
E c <C uj <§; t _1 is finite and approximately frequency- 
independent. In view of the close connection between 
the non-linear conductivity and the Fock contribution to 
the equilibrium persistent current discussed above, we 
think that our result is physically more reasonable. The 
fact that for \lo\ 3> E c the non- linear conductivity (^5[) is 
exponentially small is closely related to the exponential 
suppression of the contribution from Matsubara frequen- 
cies larger than E c to the average Fock current Jf in 
Eq.©. 



C. Averaging over disorder 



IV. POSSIBLE EXPERIMENTAL TESTS 



According to Eqs.(p2|) and ((53|) the average non-linear 
conductivity can be expressed in terms of the flux- 
dependent part of the average linear conductivity. The 



For simplicity, let us consider a time-dependent but 
spatially constant external electric field along the circum- 
ference of the ring, 



G 



E(t) = E(lu) cos(wi) 



(59) 



Because at zero wave-vector and finite frequencies the 
polarization 11(0, oj) vanishes, this field is not screened. 
Hence, q • E(q = 0,u/) in Eq.p7|) can be identified with 
the external field E(oj). Experimentally, the field ( |5g ) 
can be generated by a time-dependent magnetic flux 
through the center of the ring, 



(p(t) = 4> + 4>(<jj) sin(u)t) 



(60) 



By Faraday's law of induction, the relation between E(u>) 
and d>(oj) is 



eLE{oj) = 2iru 



0(w) 



(61) 



Note that in the experiment g the current was mea- 
sured in the presence of such a time-dependent flux 
with frequencies in the range between 10 and 10 3 Hz. 
In this range no frequency-dependence of the current 
was detected, so that apparently the measurements were 
performed in the static limit. Note, however, that in 
principle one should distinguish between the thermody- 
namic equilibrium current that is determined by the flux- 
dependent part of the free energy, and the dynamic cur- 
rent that is obtained from the time-dependent response in 
the limit of vanishing frequency [^lj . In the present work 
we are interested in the frequency range A < w < t _1 , 
corresponding to frequencies between 10 s and 10 13 Hz. 
We predict that the static non-equilibrium current should 
become exponentially small as soon as the frequency of 
the electric field exceeds the Thouless energy E c . This 
effect can be used to directly measure the Thouless en- 
ergy of a mesoscopic ring. For the rings used in Ref. 
the time-dependent non-equilibrium current should dis- 
appear for ui « 10 10 Hz. We would like to emphasize that 
this prediction can be verified without any modifications 
of the experimental setup used in Refs. [^J^]. 

Because the non-equilibrium current / is driven by 
an external time-dependent flux </>(w), it can be easily 
distinguished from the thermodynamic equilibrium cur- 
rent. Let us now discuss the expected size of this non- 
equilibrium current. Given the fact that the external field 
in Eq. (|5S|) has a cos-dependence, the experimentally mea- 
sured non-equilibrium current is determined by the real 
part of the non-linear conductivity, 



Re 



OJ 



2tt0(w) 



-i 2 



(62) 



For frequencies \oj\ ;$ E c we may use Eq.(|5q) to simplify 
the non-linear conductivity, so that in this regime we 
obtain after some rescalings 



f( 2 ) 



cE c 

4>o 



2tt</>(w) 



l^r 1/2 /(0,^,T) , (63) 



where <t> = 4>I4>0, u = u/E c , 7 = T> / (L^\lj\) . Defining 
the dimcnsionlcss variable 



x = 2 sin(27T(ft) 

the function /(</>, tu,"f) can be written as 

Xfr + X 2 } 



f W^n) = 8cos(2tt<; 



1 + b + x^y 



(64) 



(65) 



A graph of /((/>, w, 7) as function of (f> = 4>/4>q for oj = 0.1 
and 7 = 1 is shown in FigM Obviously the size of the 

_(2) _ 1—. _ 

current / in Eq. (|63|) is determined by three experi- 
mentally controllable parameters: 4>{oj) , (j), and oj. Let 
us find the values of these parameters that maximize the 
current. Obviously 4>(oj) should be chosen as large as pos- 
sible. It should be kept in mind, however, that Eq. (p3[) is 
the quadratic order in a systematic expansion in powers 
of 4>(oj) j 4>o ■ Higher orders should be negligible as long 
as \<j)(uj)\ <C 0o- Thus, the largest value of \4>(oj)\ where 
Eq. (p3|) can be expected to be accurate is 



Z7T 



(66) 



Next, consider the optimal choice of the frequency. Be- 



cause of the factor of 



1-1/2 



in Eq.(|63|), it is advanta- 



geous to choose the frequency as small as possible. How- 
ever, our perturbative calculation breaks down for fre- 
quencies of the order of the mean level spacing A. The 
optimal choice is therefore 

u w A . (67) 



Finally, from Eq.(65) and Figs.|5|, || it is clear that the 
function /(X, 7) = //cos(27T0) has a maximum at a 
value X m (i) = O(l) (provided 7 is not too large). This 
implies that for \oj\ C £ c the non-equilibrium current is 
maximal at the static flux <fi = <fi m , where 



(68) 



The function X m (j) is shown in Fig.g. From Eq. (J64|) 
we also see that the width A</> of the interval around 
4> m (modulo <po/2) where the current is enhanced is 
A(j) w (/) |(I)| 1 / 2 /(47r). Outside this interval, which is 
rather narrow for \uj\ <C 1, the non-equilibrium current 
is much smaller than at the maxima (see FigM. In fact, 
for — <p m \ A(j) the parameter \X\ in Eq.(J64|) is large 
compared with unity, so that we may approximate 




7) 



|w| 5 / 2 cos(2tt0) 
4 sin 5 (27T0) 



\X\ » 1 , 



(69) 



where we have assumed that 7 ;$ 1. Using 2 sin 2 (2) 
1 — cos(2x)], the non-equilibrium current in this regime 
can be written as 



cE c 

4>o 



2tt</>(w) 



sin(47T(/)) 



[1 - cos(47r</>)] 3 



(70) 
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For the parameters of the experiment M (taking now 
the spin degeneracy into account), we find that at the 
optimal values of the parameters given in Eqs.(|6(] |68| ) 
the maximal amplitude of the non-equilibrium current 

is I max ~ 5 x lO^evp / L. This current has the same 
order of magnitude as the equilibrium current measured 
in Refs. [^],[|, an d therefore should be measurable with 
the available technology. The rather pronounced peaks 

—(2) 

of I as function of <f> for small values of (j) (modulo 
</>o/2) distinguish the non-equilibrium current from the 
thermodynamic persistent current. 



external flux <f>(t) of the form (|6(j) , with frequencies in the 
range A ,$ w < t~ . We hope that such an experiment 
will be done in the near future. 
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V. SUMMARY 



In this paper we have presented a thorough theoreti- 
cal analysis of the static non-equilibrium current 1^ in a 
mesoscopic metal ring threaded by a magnetic flux, which 
is generated in quadratic response to a time-dependent 
electric field. Using a path integral approach, we have 
shown that this non-equilibrium current is closely related 
to the screened exchange contribution to the thermody- 
namic persistent current |l3| , |i^| . In fact, from Eqs.(^9|), 
(p|), and (|||) it is obvious that the weight of the RPA 
screened interaction in the exchange correction to the 
equilibrium persistent current at fixed energy- momentum 
transfer is essentially given by the non-linear conductiv- 
ity er( 2 ) (ui, 4>) associated with I^ 2 \ This observation has 
allowed us to derive the relation ( ^8| ) between the non- 
linear conductivity and the flux-dependent part of the 
polarization. Thus, the disorder average a^'iu)^ <f>) of the 
non-linear conductivity is directly related to the weak- 
localization correction to the average polarization, which 
in turn can be expressed in terms of the weak- localization 
correction to the linear conductivity. In other words, the 
average non-linear conductivity is directly related to the 
weak localization correction to the linear conductivity. 

A measurement of the average static non-equilibrium 

_(2) 

current I in a mesoscopic metal ring as function of 
frequency and static flux would be very interesting for 
several reasons. Such a measurement would directly 
probe the weak localization corrections to the frequency- 
dependent polarization and the linear conductivity of the 
system. In contrast to Kravtsov and Yudson [Q , we pre- 
dict that the static non-equilibrium current should be ex- 
ponentially suppressed when the external frequency ex- 
ceeds the Thouless energy. This is in agreement with 
the close connection between the non-equilibrium current 
/( 2 ) and the Fock contribution Ip to the equilibrium per- 
sistent current. The latter is known to be exponentially 
suppressed if the temperature becomes larger than the 
Thouless energy fL3f . Note that in this case the tem- 
perature acts as an infrared cutoff, just like the external 
frequency in the case of the non-equilibrium current. Our 
theory can be verified experimentally from the measure- 
ment of the time-independent non-linear current response 
of a mesoscopic metal ring pierced by a time- de-pendent 
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FIG. 1. (a) Graphical representation of the cur- 

rent-density jq in linear response, see Eq.(^). The solid 
lines represent non-interacting Green's functions for fixed dis- 
order potential, the wavy lines represent the fields V q , and the 
black triangle denotes the current vertex (k x + a + q x /2)/m. 

(2) 

(b) Graphical representation of j q . 





FIG. 2. Hartree contributions to the functional average of 
the non-equilibrium current density given in Eqs.(^3|) and 
(p5|). (a) Contribution from the linear response current 
(b) Contribution from the quadratic response j q 2 ^ . The dou- 
ble wavy line is the RPA interaction. 




FIG. 3. Fock contribution to functional average of j^ 2 \ 
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FIG. 4. f(<j},ui,-y) as function of 
7=1. 



o for lj = 0.1 and 
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FIG. 5. f(X, 7) = //cos(27T0) as function of X (see 
Eq.(||)) for 7 = 1. 
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